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The temperature dependence of the order parameter of the Schwinger model is 

calculated in the euclidean functional integral approach. For that we solve the 

model on a finite torus and let the spatial extension tend to infinity at the end of the 

computations. The induced actions, fermionic zero-modes, relevant Green functions 

and Wilson loop correlators on the torus are derived. We find the analytic form of 

the chiral condensate for any temperature and in particular show that it behaves 

like {tyfy) ~ — 2Texp(— ir^/nTfe) for temperatures large compared to the induced 

photon mass. 
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1 Introduction 

The study of exactly soluble field theories has always received a good deal of attention in 
the hope that they might shed some light on more realistic theories. One such model is 
the Schwinger model or quantum electrodynamics for massless fermions in two space-time 
dimensions [1]. It serves as an important tool in illustrating various (related) field theo- 
retical concepts such as mass generation, dynamical symmetry breaking, charge shielding 
or fermion trapping [2]. 

More recently, the Schwinger model and its extensions have been used as toy models for 
studying Baryon number violating processes [3]. Since the semiclassical expansion in the 
standard model breaks down at the interesting energies or temperatures, it is important to 
compare the realistic calculations to the corresponding ones for simpler models. As simple 
soluble models which, for the problem at hand, have features identical to the standard 
model, the non-linear sigma model [4] and the Schwinger model (coupled to a scalar field) 
[3] have been considered. 

The influence of an external gravitational field on the ultra-violet and infra-red prop- 
erties of the model has also been investigated [5]. This problem is not completely solved 
up to now but the partial results indicate that qualitatively all long-distance features of 
the flat-space-time solution persist in curved backgrounds. 



The zero-temperature Schwinger model has been solved some time ago by using oper- 
ator methods [6] and more recently in the path integral formulation [7]. Some properties 
of the model (e.g. the non-trivial vacuum structure) are more transparent in the operator 
approach and others (e.g. the role of the chiral anomaly) are better seen in the path 
integral approach. 

When studying the finite temperature 0(3)-sigma model coupled to fermions (as a 
relevant model for the B- number violating processes in the electro- weak theory [4]) we 
realized that no satisfactory solution of the simpler finite-temperature Schwinger model 
is available, at least in the path integral approach. Therefore we decided to investigate 
the Schwinger model at finite temperature and in particular determine the temperature 
dependence of the chiral condensate by using functional integral methods. To study the 
temperature effects we solve the Schwinger model on the two-dimensional torus to avoid 
the infra-red divergences, which sometimes plague the analysis of two-dimensional gauge 
theories. 

Previous finite- volume solutions of QED 2 have clarified various aspects of the model. 
For example the role of the nonintegrable phase factor exp(ie J A^dx^) has been under- 
stood by studying the model on a cylinder, that is by considering the zero-temperature 
model where the spatial dimension is given by a circle [8]. Recently the Schwinger model 
has been solved in the path integral approach on the 2-dimensional sphere and the role 
of the fermionic zero modes has been emphazise [9]. Also it has been considered on a 
2-dimensional disk assuming the most general (local) self-adjoint boundary conditions for 
the fermions [10]. 

However, these geometries are not the proper ones to study finite temperature effects. 
For that we must assume that the fields are ( ant i) -periodic in the time direction with 
period /3 = 1/T. In order to justify the manipulations below (in particular the treatment 
of the Grassmannian integrals) we prefer to work on a finite space-time volume such that 
the relevant wave-operators possess discrete spectra. Thus we assume that the spatial 
extension of space-time is finite as well and consider the Schwinger model on a torus 
[0, f3] x [0, L] with volume V = (3 ■ L. We assume periodic boundary conditions in the 
spatial direction. Only at the end of our calculations do we let the spatial extension L 
tend to infinity. 

The only other solution of the Schwinger model on the torus we are aware of, besides 
assuming periodic boundary conditions for Dirac-Kahler fermions instead of finite tem- 
perature ones, concentrated mainly on finding and discussing the fermionic zero-modes 
[11]. Also, our trivialization of the [/(l)-bundle over the torus is different from the one in 
[11] but is very convenient for determining the explicit spectrum of the Dirac operator. 

Our main results are the analytic form of the (zero-mode truncated) effective actions 
and zero-modes in the various topological sectors, the two-point Green function of A 2 — 
m 2 A on the torus and as a consequence the explicit formulae (5.9), (5.10) and (5.12) 
for the temperature dependence of the chiral condensate (\I/\I/). We find that (fyty)? is 
equal to — (m 7 /27r) • e 7 for T = 0, in agreement with the zero-temperature solutions 



[7,9], and tends to — 2Texp(— 7rT/m 7 ) for temperature large compared to the induced 
photon mass m 7 = e/y/n . Our high temperature result differs by a factor 2 from the 
only other (formal) finite temperature path integral solution we are aware of [13]. Also, 
the torus calculation nicely splits the chiral anomaly responsible for the non-vanishing 
condensate, into its high- and low energy parts. The first (coming from the gauge invariant 
ultraviolet regularisation) only affects the numerical value of (\E , \[ / ) , whereas the latter 
(coming from the chiral asymetry of the zero-energy sector) is responsible for the non- 
zero fifty). This is the reason why ordinary perturbation theory, although leading to 
the correct f/V-contribution via the wellknown anomaly graph, fails to predict a non- 
vanishing condensate. At the end we obtain analytic expressions for (thermal) Wilson 
loop correlators and in particular their temperature dependence as well as the gauge 
invariant fermionic two-point functions. 

2 The Schwinger Model 

In this section we recapitulate briefly the euklidean Schwinger model [7], using functional 
methods. We emphasize the less known features involved when considering the model on 
a torus instead on the plane. The Schwinger model is defined by the action 



S[A,V,V] = ^ j d 2 xE 2 - J d 2 x^ipty, (1.1) 



where 

E = F i = d Q A x - 8 X A Q and p = 7 "(d„ - ieAj . 

Note that in two space-time dimensions the electric charge has the dimension of an inverse 
length, contrary to the situation in four dimensions. The model is superrenormalizable 
and requires no infinite renormalization besides a trivial redefinition of the zero energy 
density. The generating functional for the Green functions is 

Z[J,r),fj] = C- 1 fvAV^Vty e -S[At,f]-s 3/ [i]+/^,(i ( ,J" + ^ + ^) ) (L2) 

where S'j/fA] contains the terms due to gauge fixing and J M resp. (fj, rf) are c- number and 
Grassman-valued currents respectively. The normalization constant C is chosen such that 
Z[0, 0, 0] = 1. First we evaluate the fermionic path integral 

Z[A;r),fj] = Iv^Vm e /^*+/^(^W) ? ( L3 ) 

and treat the 'photon'-field as external field in this first step. 

On the torus the Dirac operator may possess normalizable zero- modes, pip = 0, and 
thus we must allow for such fermionic zero modes in the path integral (1.3). Let us assume 
that there are \k\ (our notation will become clear later on) normalizable zero-modes ip p . 
In addition there are an infinite number of excited modes ip q , q > \k\ with eigenvalues 
\ q ^ 0. The number of zero-modes can be determined as follows: 



First we note that since 75 anti-commutes with the Dirac operator the excited modes 
come in pairs 

iplpq = Xql^q =► %p {l^ q ) = -\{l^q), (1-4) 

so that ip q and 'j^ipq have eigenvalues X q and — X q and are orthogonal to each other. 
It follows that the contribution of the excited states to the (supersymmetric) partition 

function 

|fc| 00 

tr 75 e'^ 2 = ^(^,75^ P ) + J2 e "*^(^«.75^)- 
1 |fe|+i 

vanishes. Thus the partition function is just the index n + — n_ = k, where n + is the 

number of right handed (chirality 75 = 1) and n_ the number of left handed (75 = — 1) 

zero modes. On the other hand, inserting the small- t expansion (see e.g. [14]) 

(xleW'lx) ~ ^ (1 + l,Et + 0(t 2 )) (1.5) 

one immediately obtains the index theorem (see e.g. [15]) 

k = — d 2 xE = — $ (1.6) 

Z7T J ZTT 

which relates the index k of the Dirac operator to the flux of the electric field. This formula 

already shows that there are zero-modes for non-vanishing fluxes $. Also note that on 

the torus the flux is quantized in integer muliples of 2n. This is really a consequence of 

the single valuedness of the fermionic wave function (cocycle condition). 

The pairing property (1.4) does not hold for the zero- modes. Indeed, they can be chosen to 

have fixed chirality since 75 commutes with the Dirac operator on the subspace spanned 

by the zero modes. Also (1.4) would not be valid for the excited modes if we would 

impose chirality- violating boundary conditions as in [10]. But (anti) periodic boundary 

conditions are compatible with the transformation ip — > 75^ and the identity (1.4) holds 

in the present situation. 

To determine n+ and n_ separately we decompose the gauge potential as 

A lt = A li -e ltu d v 4> where = —(E - -$) + c =► E = - + A0, (1.7) 

that is into a global "instanton"-type potential A with constant field strength, E = $/V 
and a local fluctuation 5A^ = —e^ v d v (p about the instanton. The subtraction of the 
constant term in (1.7) is necessary since on the torus the Laplacian has the constant zero- 
mode and is invertible only on functions which integrate to zero. Using 7^75 = —ie fiu , ~f u 
it is now easy to see that 

which shows that the number of fermionic zero-modes is independent of <fi and hence is 
the same for A and A. But for the instanton potentials A 

p 2 = D 2 + l5 ^ (1.9) 



~ 2 

and since —p is a non-negative operator, all zero modes are either right- or left-handed 
for non- vanishing fluxes. Only in the zero-instanton sector $ = can it happen that there 
are both right- and left-handed zero modes. For example, for a vanishing gauge potential 
there is one right-handed and one left-handed zero mode. 

What we have shown then is that for $ 7^ there are exactly \k\ — |<&|/27r zero modes. If 
the flux is positive they are all right-handed, else they are all left-handed. 

After having counted the number of zero modes we proceed by expanding the 'electron'- 
field in an adapted orthonormal eigen-base as 

|fc| 00 

1 lfel+1 

and similarly \I/, so that 

(V, *) = ^2% l/> p )Op + ^2(fj, ^q)Pq (*, V) = Yl ®P$P> V) + Y1 Pqi^q, V) 

split into a zero-mode part and an excited part. Inserting this decomposition into (1.3) 
and using D^fD^f = T>aT> aDflDfi the Grassmannian integral over the a's can easily be 
done since the action does not depend on them. This way one finds for the zero-mode 
contribution to (1.3) 



VaVa 



e (.T),ip p )a p e a p (ipp,ri) 



(r},ip p )(ip p ,r)). (1.10a) 



The remaining /3-integration is performed by shifting the /3's (and similarly the /3's) ac- 
cording to 

After this shift the /3 integration yields 

/oo 
VBVB TT e X oPipQ +< < f i, , >Pq)Pq+Pq{$q,V) — g~ I v(x)G e (A;x,y)r)(y) _ feft! Up\ (1.106) 

|fc|+l 

where det' is the determinant with the zero-eigenvalues omitted and 



:x 



Ge( ^) = £«M. (l.lla) 

|fc|+l q 

is the Green function on the space orthogonal to the zero modes. Clearly this function 
obeys the differential equation 

1*1 
ipG e (A;x,y) = 5(x-y)- P(x,y), where P(x,y) = ^ ^p&Wpiv) (1.H&) 



is the projection kernel on the zero-mode subspace. Due to the pairing property of the 
excited modes the 'excited' green function anti-commutes with 75: 

7 5 G e (A;a;,i/)75 = -G e (A;x,y). (1.11c) 

Inserting now (1.10) into the path integral for the (fermionic) partition function (1.3) we 
end up with 



Z[A;r},r}] 



(f},if> p )(il> p ,ri) e~ * iW G *( AiX > y W») det' (ip). (1.12) 



Thus to determine Z(A; fj, rj) for a fixed potential we need to compute the \k\ zero- modes of 
the Dirac operator, the determinant of ip with the zero-modes omitted and the 'excited' 
Green function. The (unnormalized) fermionic 2n-point functions for a given A-field are 
now obtained by differentiation with respect to the external currents 

5 2n 
= ^(y„)^ a -(x n )...(J^(yi)(J^(x 1 ) Z[i4; ^ 77]l ' J=,F=0 (L13) 

Using (1.12) in this formula we can immediately read off that 

1. The fermionic partition function 

Z[A;0,0] = det(ip) (1.14) 

which enters in the normalization constant C in (1.2) is only non-zero if the Dirac 
operator possesses no zero modes that is for gauge field with vanishing flux. 

2. The 2-point functions are non-zero only if A admits either no or one zero-mode, 
that is for $ = or $ = ±2tt and then 

[mMer*»*9 a (zmv) = h G ^ X ^ det ^ $ = ° (1.15) 

J \-^t(x)^(y)def(»f) |*|=2rr. 

From (1.11c) it follows that the expectation values of i &P±ty, where P± = 1(1 ±75) 
are the chiral projections, are non-vanishing only for $ = ±2tt in which cases the 
Dirac operator has one zero mode ip of chirality ±1, and then 

/ V^Vm e /*^* ^>(x)P±m(x) = -tr^^P^ix)) det' (ip). (1.16) 

The higher 2n-point functions are obtained similarly. For example, for the 4-point 
functions only gauge potentials with no, one or two zero modes contribute. In 
particular for the expectation value of the operator with (chiral) charge zero 

O(tf) = (^P + ^)(^P_^), 0(e° 75 ^) = 0(f) (1.17) 



only the trivial sector with no zero- modes yields a non- vanishing result. This ex- 
plains why earlier solutions of the Schwinger model on the plane (in which the 
non-trivial sectors have been neglected) yield the correct result for this expectation 
value and thus (via clustering) the correct result for the chiral condensate. 

Note that the excited fermionic Green function does not appear in the expectation values 
(1.16), and this simplifies the computation of the chiral condensate considerably. More 
generally, the expectation values of the charge ±ra operators (^fP±^f) n are non-zero only 
for potentials with fluxes $ = ±27m, and then only the zero modes and effective actions, 
but not the excited Green functions occur in the expectation values. 

3 Effective Actions - Local Part 

In this section we compute the ^-dependence of the effective actions det' ip for the general 
gauge potentials (1.7) by integrating the chiral anomaly. Since there are some subtleties 
when zero-modes are present we derive the local ^-dependent part of the effective actions, 
although some of our arguments are not new [16]. 

To determine the ^-dependence of the effective actions we consider the one-parametric 
family of Dirac-operators 

p a = e^ a pe^ a ^p a = { 7 50, Pa} C= ^) (2.1) 

which interpolates between p and p , and calculate the variation of the zeta-function 
regularized determinants [17] 

1 1 A 

logdet'(i# a ) = -logdet'(-f I) = --— C^ 2 (s)| s =o. (2.2a) 

The zeta-function is defined as as the analytic continuation of 

oo 

CW*)=Z>7'» R(*)>1, (2.26) 

|fc|+i 

where fj, q = A^ are the a-dependent excited eigenvalues of the squared Dirac operator, to 
the whole complex s-plane. Using (2.1) in the Feynman-Hellman formula 

X q = (ip q ,ipip q ) = 2A g (Vv75<M?) 
the variation of the zeta-function can be written as 

, oo 

-fc<ipi ( s ) = ~ s Yl ac -1 /^ = - 4s $^rwv^5<M?)- 
ifci+i 



This can be further rewritten as a Mellin transform 

a ^ J |fc|+i ^ s ' J 

where the trace is only to be taken over the excited states. Inserting the asymptotic 
expansion (1.5) from which we must subtract the projection density P a (x,x) = P a {x) 
(see (1.11b) onto the zero modes we find 

^C^ (s) U=o = ~ J d 2 x Ec/> + 4 J tr (P a {x) lb <f>{x)) . (2.3) 

To integrate over a we observe that the zero-modes of p a and p are related as 

4 a) = e^^ p , p=l,...,\k\ (2.4) 

with the negative sign in the exponent for right-handed and the positive sign for left- 
handed zero-modes. These modes are in general not orthonormal, and the normmatrix 

A/- pr (a) = «\^ a) ) (2-5) 

is not the identity so that the projection density reads 

pr 

Now it follows from (2.4) and (2.5) that 



_d_ 
da 



log det (A/») = -2 I d 2 x tr {P a ^(f) (2.6) 



and this formula can be used to integrate the anomaly equation (2.3) from a, — to a, — 1. 
Together with (2.2) we end up with 

det' {ip) = det -^ det' (ip) exp (— / E<j>) (2.7) 

M x 2n J ' 

which is the (almost) factorization of the determinants into a global, A-dependent, and 
local, ^-dependent, part. The factorization is not complete since det N = detjV(l) still 
contains a coupling between the instanton potential A and the local fluctuation 5 A via 
(2.4). Since E in (1.7) depends on both $ and 0, there is also an apparent coupling in the 
last factor in (2.7). However we shall see later that not all functions (f> are permitted and 
that for the allowed ones this factor does not depend on $. Also note that the local part 
does not depend on the geometry of the torus. In particular the induced photon mass 
m 7 = e/ ' \pii is the same as on the infinite plane. 



4 Effective Actions - Global Part 

With the factorizations (2.4) and (2.7) of the zero modes and effective actions into local 
and global factors the problem of finding the 2-point functions (1.16) (or 2n-point func- 
tions of charge ±ra operators) reduces to the problem of computing the zero modes and 
effective actions for the instanton potentials A. For that we need to know the explicit 
form of these potentials. 
An instanton potential can always be decomposed as 

A\ = -— x 1 + -^-h + d \ and A 1 = -^-h 1 + d 1 X (3.1) 

with constant h^. The least trivial term proportional to $ yields the constant field 
strength and corresponds to a particular trivialization of the £7(l)-bundle over the torus. 
In other words, the gauge potentials at (rc^x 1 ) and (x°,x 1 +L) are necessarily related by 
a non-trivial gauge transformation 

A^x , x 1 + L) — ^(a; ,^ 1 ) = d^a, where a = — — x°. (3.2a) 

P 

Accordingly the fermionic wave functions transform as 

^(x\x l + L) = e'XxV 1 ). (3.26) 

These boundary conditions must be supplemented by the finite temperature boundary 
conditions in the x°-direction 

A li (x°+p,x 1 )=A li (x°,x 1 ) and tp(x° + (3,x 1 ) = -ip(x ,x r ). (3.2c) 

The constant terms in (3.1) are the harmonic pieces of the gauge potential. Even in the 
trivial sector $ = they yield the (dynamical) nonintegrable phase factors 



ifA^dxV 



exp \2Tri(h no + hini)] (3.3) 



for loops which wind no-times around the torus in the x° direction and ni-times in the 

x 1 -direction. 

The last term in (3.1) is a pure gauge term and will drop out at the end of the calculations. 

Now we consider the sector with $ = and those with $ ^ in turn: 

4.1 Trivial sector $ = 

~ 2 

In this sector —If) = (V—h) 2 Id is simple and possesses the double degenerate eigenvalues 
(of course, the spectrum is not affected by the gauge term in (3.1)) 

27T 2 27T 2 1 

fJ-m = (-o-) (mo - a ) 2 + (— ) (mi - a x ) 2 where (a , a x ) = (- + h , hi). (3.4) 



The corresponding zeta-functions £(s) = Y^ fi m s are the generalized Epstein functions [18]. 
Applying the Poisson-resummation formula the derivative at s = can be expressed in 
terms of Jacobi theta-functions as [19] 



d 
ds 



C(*)l 



s=0 



-2 log 



77(ir 



-e 



2 + a 
2 -Ol J 



(0,*r) 



where r = L//3, 



is the ratio of the two circumferences of the torus and 



ri(iT)=q 1/M Yl(l-q r - 



q = e 



-2lTT 



(3.5a) 



(3.56) 



n>0 



Dedekind's eta-function. We have adopted the conventions in [20] for the theta-functions: 



e 



a 
lb 



(Z,IT 



E' 



,— iTT(n+a) 2 +2iri(n-\-a)(z+b) 



(3.5c) 



Taking the degeneracy of the eigenvalues into account, the zeta-function regularized de- 
terminants (2.2a) in the sector with no zero modes are 



det(ip) 



rj(iT 



-e 



2 + a 

- — a\i 



(0,zr) 



(3.6) 



4.2 The non-trivial sectors 

We have seen that the excited eigenmodes of the Dirac operator come in pairs with 
opposite eigenvalues. Since 75 commutes with the squared Dirac operator the chiral 
projections P±ip q of these modes are eigenmodes of — Jj) 2 . Thus the excited eigenmodes of 
the squared Dirac operator come in pairs as well and two partners have the same energies 
\x q = \ 2 but opposite chiralities. We have also shown that there are exactly \k\ — |$|/27r 

chiral zero modes of the Dirac operator and hence of —]p 2 . Since Ip in (1.19) differs only 

by the constant 2<&/V in the two chiral sectors, the \k\ zero modes with chirality, say +1, 

are at the same time excited modes with energies 2$/V and chirality — 1. Due to pairing 

~ 2 
there are also \k\ excited modes with chirality +1, and so on. Thus —Tp possesses the 

following spectrum 

degeneracy = |$|/27r 

2n|$|/V degeneracy = |$|/7T 

Note that, contrary to the situation in the zero-instanton sector, the spectrum does not 

depend on the harmonic part h of the gauge potential A. 

The zeta-function is now proportional to the ordinary Riemann zeta-function and one 

obtains 

/7rV\ I*I/4tt 

det'(^) = (— ) (3.8) 



Hr 



for the determinant of the non-zero eigenvalues. 



10 



Inserting (3.6) and (3.8) into (2.7) we end up with the following formulae for the 
determinants in the sectors with \k\ zero modes: 



det(ip) 



r}{ir\ 



■e 



exp (^ / 0A0) k = 

det'(^)=(^) |fe|/2 det^ exp(^|0A0+^|0) k* 0, 

where is the theta-function in (3.6) and we have expressed -E in terms of <fi and $ (see 
(1.7)) and made the dependence on the dimensionful electric charge explicit (recall that 
m 7 = e/yfjr). 

5 Computing the Zero-Modes 



Z. 2 



According to (2.4) we only need to evaluate the zero-modes of —If) . Also we shall only 
consider expectation values of gauge invariant operators and thus may set the gauge part 
in (3.1) to zero. Since the gauge potentials (3.1) do not depend on time the eigenmodes 
should be proportional to exp(ipoX°). The finite temperature boundary conditions (3.2c) 
then require that the momentum is quantized as po = (2p — l)ir/(3 with integer p. If we 
further eliminate the /ii-part of the instanton potential we are lead to the following ansatz 

= i(2p-l)7nr°//3 i2-wh 1 x 1 /L 



% = e *(2p-l)**70 e i2nh lX "/L £ / h 



Z 2 



for the zero modes. Inserting this ansatz into the zero mode equation Jp \ v = yields 
,d 2 $ 2 2 |$k, n , 1 L. 

i^f-^y +-r7-)^ = ' where y = x +^(p- a o) 



which shows that £ p is the ground state wave function of a harmonic oscillator and thus 

£ P = exp 



M^ + ICp-oo)} 2 !, (4.1) 



2V l k 

where we have used the index theorem (1.6). These functions do not obey the boundary 
condition (3.2b) but the correct eigenmodes can be constructed as superpositions of them. 
For that one notes that 

Xpi^y+L) = e-*** /* e**i xp+kix^x 1 ) (4.2) 

so that the sums 

& = (^) llA E ^ nhl X P+nk , P = 1, .... |*| (4.3) 



11 



obey the boundary conditions and thus are the \k\ required zero- modes of the Dirac 
operator. The overall factor normalizes these functions to one (the norm can be computed 
by using 4.2)). Also, modes with different p are orthogonal to each other, so that the 
system (4.3) forms an orthonormal basis of the zero-mode subspace. These functions can 
be written in terms of Jacobi theta functions as (we drop a constant phase) 



i> P 



( 2\k[ 

{ fi 2 V 



\ 1 / A ^,m[h x /fi-kx x 1 /V\ q 



x 1 jL + (p-a )/k 
kx°/(3 + oi 



(0,i\k\r) 



(4.4) 



where r has been introduced in (3.5a). Since these modes are orthonormal the determinant 
of the norm matrix M in (3.9) is just one. The unnormalized zero modes which enter the 
normmatrix M can now be obtained from the zero modes (4.4) by the transformation (2.4) 
with a — 1. These completes our discussion of the zero-mode sector. We have computed 
all \k\ zero modes of If) explicitly in terms of Jacob i-theta functions. 
For the two-point functions of interest (1.16) only the sector with one zero mode con- 
tributes, and since then the normalized zero mode is just ip/v^f, the normalization con- 
stant cancels with M in (3.9) so that 



WW ef* ip * V(x)P+m 



1 



x) 



-i\®±\ e 



2 =F2</> 



cxp 



m: 



where 0± are the theta-functions in (4.4) with p = and k = ±1: 



+ 



ek 
V 



, (4.5a) 



0+ = 



x 1 J ' L =p oo) 
x°/p ± Oi 



(0,ir) 



(4.56) 



(we used that the modulus of this function is unchanged if the sign of the " lower param- 
eter" is changed) and the ± sign refers to the chirality of the zero modes. 
Our trivialization of the U(l) bundle which leads to the form (3.1) of the instanton 
potential differs from the one chosen in [11] and so do our zero modes. But one can show 
[12] that the modes in [11] are obtained from our modes (4.4) by the corresponding change 
of trivialization, as required. 

5.1 Calculating the bosonic path integral 

After having solved the fermionic integration we are now left with the functional integrals 
over the A-field. In particular we shall evaluate the two point functions 



(VP±9) 



fVAe-y E2 fWW e^ i ^m{x)P ± m\ 



x) 



JVAe-lf E ' 2 JV^Wef^* 



(5.1) 



where E is given in terms of and $ in (1.7). Clearly, since the integrands are expressed 
in terms of {0, /i M , $} it is natural to change the integration variables from A to these 
variables. 
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First one notes that there is a one to one correspondence between E and {A0, $}. This 
becomes a one to one correspondence to {0, $} if we demand that is orthogonal to the 
kernel of A, that is if it integrates to zero. Second, for given {0, $} the nonintegrable 
phase factors (3.3) are in one to one correspondence to the phases exp(2i7ih^), and thus 
to the h^ modulo 1. Furthermore, from dA = AA we see that there is also a one to one 
correspondence between the divergence of A and functions A which integrate to zero. To 
summarize: we have shown that the transformation 

A n — ► {& A, /i„, $} where (j, = \ = and < h„ < 1, (5.2) 

defined in (1.7) and (3.1), is one to one. 

Next we need to calculate the Jacobian of this transformation. This is conveniently done 

by expanding all fields in eigenmodes of —A, e.g. 



(f)(x) = ^—^<j) m e 2«r(™o*7P+mi*7^ > (5.3) 

where m = (mo, mi) = is excluded because of the constraint on in (5.2). In terms of 
the coefficients the transformations (1.7) and (3.1) read for $ = 



A 0: m\ = 2i7r /-mi rm \ /<j> n 
Ai m ) L \rm m ± J \X n 



(m^O) 



Ao '° i = 2tt i ^ ° i r° 

A h0 ) n I 1/y/Tj Ui 



and the Jacobian of this transformation is just J = (27r) 2 det'(— A) and thus independent 
of the dynamical fields. In the nontrivial sectors labelled by k we can write A^ = — $/V • 
x 1 S^ 1 q + 5A^ and the above transformation applies then to SA, so that 

i 

/ va^ = j2 v8A » = j Yi I dh ° dhi I n d ^ dx ^- ( 5 - 4 ) 

Since in expectation values the same Jacobian appears in the numerator and denominator 
we may neglect J in what follows. Also, gauge invariant operators (like in (5.1)) are 
independent of the gauge function A and the integration over the A m drops in expectation 
values as well (in our decomposition A = corresponds to the ghost free Lorentz gauge). 
Inserting (3.9) and (4.5a) into (5.1) we find 

( »i» ± ») = zl e -^v K/>|e ± iy^e-^*> 



where 



r[0] = ^y"0(A 2 -m2A)0 (5.56) 
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and one integrates over which integrate to zero. Q± and O are the theta functions in 
(4.5b) and (3.6), repectively. Note that the last term in the exponent in (4.5a) is absent 
for the allowed 0's. 
Finally, using 

<fh\e±\ 2 = fd 2 h\e\ 2 = J— (5.5c) 

and performing the Gaussian functional integral yields 

(VP ± V) = JWLh e -2«*/e*V e 2e*K(x,x) where K f xx \ = l x \ 1 U\ ( 5 . 6 ) 

p A 2 — m^A 

is the Green function on the space of functions which integrate to zero. Using the eigen- 
modes (5.3) of —A and observing that the excited modes span the space of functions 
permitted in the path integral (5.5a), we obtain 



K( x ,y) = Y. *¥fi" ] =* *('■*) = 4v E (~ 1-.): 

^— ' n z 4-m z ii — m A \/ ^— ' K ii— ii— -4- m z ' 



4>\n{ x )4>m{y) ^ 

H^ + m^m -"••-< m 2y Z^K j2m /im + m 2 



where the eigenvalues fj, m of —A are the ones in (3.4) with a^ and a\ set to zero. Using 
the identity 

± ^■^'"J- 2 *!) (•>„) (5.7) 

^— ' a +mi a sinh(7ra) 

mi=-oo - 1 - v ' 

for $ = the summation over m\ can be carried out and one obtains (n = m 1 ) 

rrtKiv ,A- l coth(Lm 7 /2) t 1 ^ r coth(n7rr) coth(7rre(n)) i 

m ^ KX ' X) ~mlV 2/3m 7 + 12 + 2n ^ I n £(n) J ( j 

where £ 2 (n) = n 2 + (f3m 1 /27r) 2 . The first sum on the right hand side can be expressed in 
terms of Dedekinds function as 

£ -(coth(vrrn) - l) = 2 £ !_ = -2 £ log (l - e "*-) = -2(^ + log^r)) 

n>0 n,r>0 r 

and thus cancels against the third term on the right hand side in (5.8) and the \r] 2 \ in 
(5.6). Also note the U-dependent exponential factor in (5.6) cancels against the first term 
on the right hand side in (5.8) so that finally 



1 



B 



7T , ,1 



qp ± q) = -- exp - coth(-Lm 7 ) e b &"*> e ~^ m ^ (5.9a 



/3m 7 2 



where 



F(x ) = v (- l ) 

£^n ^^ + {x/2^ j 



H{x,r) = Y y - = ; . (5.96) 

^ Vn 2 + (X/27T)2 e Ty/(2«n)*+ X * _ 1 
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The formula (5.9) is the exact form of the two point functions we have been aiming at. 
This formula simplifies considerably if we let L — > oo, in which case H(.,.) — > and 
coth(.) — > 1 so that 
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(^) = _^exp ^__ e F(/3m ^, (5.10) 

p i pm^- 

where we have used that P + + P- = Id. 

The formula (5.9a) for the finite temperature chiral condensate on a finite intervall and 
the corresponding formula (5.10) after the limit L — > oo has been taken, are the main 
results of this paper. Let us now derive the low - and high temperature limits in turn: 
To study the low temperature limit of the chiral condensate we regularize the indiviudal 
sums in (5.9b) as [21] 

V—— ~ - + 7 + 0(s) 
A^ n i+s s v ) 

_ i i °° 

V J— - - log ^ + 2 V K (lx) + 0(s) (5.11) 

where 7 = 0.57721... is Eulers constant. The difference is finite as s tends to zero and 
inserting it into (5.9) and (5.10) yields 

00 

(^) = _I^2 e 7 e 2/(/3m 7 ) where U x )= f 1 -^ 7? rdt, (5.12) 



where the last integral represents (up to a sign) just the last sum in (5.11). This altena- 
tive representation of the chiral condensate for L = oo has been derived in the bosonised 
Schwinger model on the cylinder in [8] . Since the integral / tends to zero for low temper- 
atures we find 

— 772 

(**) — > — -^ e 7 for T -> 0, (5.13) 

2n 

and this is the well known result for the zero temperature Schwinger model [7]. 

To study the high temperature limit we observe that F(x) tends to zero as x ^C 1 or if 

the temperature is big compared to the induced photon mass. In this region 

(M) — ► -2T e -* T/m ~< for T ^ oo. (5.14) 

Thus the chiral condensate vanishes exponentially fast for T ^> m 7 . 

For the intermediate temperatures < T ~ m 7 we must retreat to numerical methods to 

evaluate the infinite sum which defines the function F appearing in (5.10). The result of 

the numerical calculations are depicted in Figure 1. 

Note how the function (^/^/)(T) resembles the behaviour of an order parameter in a system 

which suffers a second order phase transition. However, the chiral condensate does not 

really vanish at any finite temperature, it tends to zero 'only' exponentially, contrary to 

the phase transition case. Thus in a strict sense the chiral symmetry remains broken at 

all finite temperature. 
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6 Further Correlation Functions 

In this section we determine the correlation functions of Wilson loops, products of the 
field strength, thermal Wilson loops and gauge invariant fermionic two-point functions. 
Clearly, for the purely bosonic operators only the zero-instanton sector yields non- vanishing 
expectation values. Since in this sector A =harmonics+e<90 one has 

<°M> = V2r jv^e-m ' ((U) 

for an operator 0[A|. Here we have used the second formula in (5.5c) and is the theta- 
function in (3.6). Hence, up to the /i-integration we remain with Gaussian integrals with 
respect to the effective action (5.5b). 

6.1 Correlation functions of the field strength 

The n-point functions of the field strength E = A<f) are particularly simple, since E does 
not depend on the harmonic pieces in the gauge potential. Because of the symmetry 
— > —(f) of T all expectation values of odd powers of E vanish. The 2n-point functions 
can be calculated from 

(E( Xl ) ■ ■ ■ E(x 2n )) = JV(f) e ~y^rM ' E{X2n) = A * • • ■ A snM*i) • ■ • <h*))- 

and applying Wick's theorem to the last expectation value allows one to express them in 
terms of the 2-point function 

(E(x)E(y)) = A x A y ((l>(x)<j)(y)) = A x A^(:r-y,m 7 ), 

where the Green function K corresponding to V has been introduced above (5.7). Since 
— A0 m = fi m (p m we easily find that 

2 

(■■■) = J2 tlWMy) [1 f^-j ] = -A.G(:r - y, m 7 ). (6.2) 

Now we use the (almost) completeness of the excited eigenmodes ra and obtain 

(E(x)E{y)) = S(x - y) - — - m*G(x-y; m 7 ) (6.3) 

where G denotes the massive Klein-Gordon propagator. This shows that E is (up to 
contact terms) just a free massive pseudo-scalar field. Using (5.7) the propagator can be 
written as 

1 ~ cosh[7rr£(n)(l-2£7L)]e 2 ™^ 1 

( ^ mi) ~ 4^ ^ e(n)sinh[7rre(n)] ^V' {QAa) 
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where £(n) has been defined below (5.8). In the limit where L tends to infinity this 
simplifies to 

1 ^ p -27rK(n)C 1 -i<°]//3 

G «^> = s £ — W) — ■ <v) 

This shows explicitely that the 2-point function of the field strength falls off exponentially. 

6.2 Wilson loops 

Since loop integrals of local Wilson loops (that is loops without windings) also do not 
depend on the harmonic pieces one finds 



J'V<f)e- m exp (ieLA(j) 

I ie§A\ V . 

where the loop encloses the region T> and we have used that the loop integral is equal to 
the flux of E = A0 through T>. Using (6.3) the resulting Gaussian integral yields 

(e ie ^> = exp ( - ^ / (E(x)E(y))d 2 xd 2 y 

V Z JVxV 

( e 2 r A 2 f i \ 

= exp ( - — A v - —^ -m 2 / G(x-y,m-,) ), (6.6a) 

V l L V JVxV J/ 

where A-p is the area enclosed by the loop. For example, for a rectangle T> one finds for 

L = oo 

(e ie ^) = exp(-^[l-m 2 | G(^;m 7 )^]). (6.66) 

For low temperatures (compared to the induced photon mass) and loops with edges 
long compared to the inverse photon mass we may use the zero-temperature propaga- 
tor G(£,m-f) = K (m 7 |£|)/27r in (6.6b) for which the integral in the exponent is, up 
to an exponentially small term, just m~ 2 . It follows then that at low temperature the 
interaction between two charges is exponentially small 

6.3 Thermal Wilson loops 

The expectation values of a string of thermal Wilson loops 

V{u) = exp (ie / A (x°,u)dx°) (6.7) 



depend on the harmonics because a thermal Wilson loop has winding number one. Since 
J A dt = 2iih Q — dif (pdx° the h- integral for a string of q loops becomes 



/ d 2 h\9\ 2 e 2 * = JX e -V/2r_ 
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The remaining (^-integration yields 

(II P(«0> = e-^ /2r exp [^ ( I d 2 xd 2 y £ 5^ - Ui ) d 2 xl K(x, y) £ o^ 1 - 



Inserting the explicit form of K (see above (5.7)) one finds the following one- and two-point 
functions 

(V(0)) = exp ( - -^ coth(-Lm 7 ) 

(P(u)P(O)) = (^(0)) 2 exp ( - ^7r/?ra 7 S(u)) (6.8a) 



where 5 denotes the function 



„ s cosh \mAL — 2\u\) 
S(u) = 2 jl ' U . (6.86) 

The normalized higher n-point functions 

nu ^ y %ff"'» ^ 

are then just products of the normalized two-point function 

V(u 1 ,...,u q ) = Y[V(u i ,u j ). (6.10) 

i<j 

Since — TlogV(ui, . . . ,u q ) is to be intepreted as the (zero energy subtracted) free energy 
of q static charges at positions u±, . . . ,u q , (6.10) shows that the free energy of q charges 
is just the sum of the free energies of the individual pairs. 
For L = oo the one and two-point functions (and thus the g-point functions) simplify to 

(P(u)) = exp ( - 7r/3m 7 /4) and V(u, 0) = exp ( - -7rf3m^e- m '< M ), 

in complete agreement with the analytically continued one and two-point functions ob- 
tained via bosonization in [8]. Note in particular that the free energy of a single charge 
is finite 

TTTfi 

AF = -T\og(V(0)) = ^ 1 (6.H) 

and that the (zero energy subtracted) free energy of two static charges falls off exponen- 
tially leading to a Yukawa force 

F = _ e - m ,,\u\ ( 612 ) 

between them. We see that the classical Coulomb force between external charges is 
shielded. 



6.4 Gauge invariant fermionic two-point functions 

The gauge invariant chiral two-point functions 

S±(x,y) = (V(x) e ie -iy A ^P ± ^(y)} (6.13a) 

may lead, via the LSZ reduction procedure, to chirality violating scattering amplitudes. 
They can be calculated in the same way as the chiral condensates (5.5a) which they must 
reproduce for coinciding points. One obtains 

s (x v) _ - 1 e -**i*v Jd 2 h4(x)My) [v± e -m-eim + m) + ieF v A 

b+{x > y) ~ ~ f<Ph\±e\*fv<j>e-m ' {b - Ub} 

where the zero mode ^q{x) has been computed in (4.4). The integration over the har- 
monics yields 

d 2 hijl(x)My) e ieJ y A = -£L e~^ x '^ e^-C^ 9 "^, (6.14) 

v 2r 



where C is some normalisation constant. The 0- integration is again Gaussian but due to 
the phase factor in (6.14) involves derivatives of the bosonic Green function: 

V<f) e -n<t>]-e(Hx)+4>(y))-iefe^d„w = c exp[K '(0,0) + K(x,y)] x 

y 
exp [_ ±jd^dr) a e;zid e d v ?K(t,r))], (6.15) 

X 

where C is again some constant. For simplicty we shall first consider the zero temperature 
and infinite volume limit only. There K depends only on the distance of its two arguments 
and the integral on the r.h.s of (6.15) vanishes identically. Hence we are then left with the 
computation of the Green function K(x,y). This Green function, which is the difference 
of the masseless- and the massive Klein-Gordon propagator (see 5.7) simplifies to 

1 —1 

K ( x > y) = ~ [ G ( r; m T = °) - G ( r ; m ~<)] = g -- ! [^sM + #o(«v)] » ( 6 - 16 ) 

where r = \x — y\ and \x is an infrared regular isat ion which would be undefined on 
the Euklidean space. Here it is fixed by demanding that S±(x,x) reproduces the chiral 
condensate (5.13). Finally substituting (6.14 -16) into (6.13) yields 

<*(*) e"/^t WI = _f^Z ^tJKptzM . (6 . 17a) 

V ^ 2tt^\x — y\ 

For large separations K^m-jr) vanishes exponentially so that the chiral two-point function 
is long range: 
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for \x — y\ approaching infinity. Note that contrary to the bosonic correlators which fall off 
exponentially due to the mass gap in the spectrum, the chiral fermionic 2-point function 
falls of like 1/ ' \fr. This is due to the zero mode fermionic states in the Hilbertspace. 

In a next step we will increase the complications by considering the case of again 
infinite volume but finite temperature. This time however we must assume the difference 
of the two arguments of the two-point function to be strictly spacelike (eg: x° — y° = 0) 
for then the integral on the r.h.s. of (6.15) vanishes again. This becomes clear if we 
notice that now only "timelike" derivatives of the bosonic Green function appear and 
these are zero for x° — y° = as can be seen using the sum representation of K(x,y) 
(above (5.6)). Hence we have again reduced the problem to the computation of K(x,y). 
However, as in the case of coinciding arguments (chiral condensate) we couldn't find 
simple expressions representing the infinite sums for finite temperature. Therefore we 
again retreat to numerical methodes to evaluate (6.15). The results are despicted in 
Figure 2. We will not reproduce here the lengthy formulas but emphasize that for non- 
zero temperature and for (x 1 — y l \ — > oo , S+(x, y) falls of like exp[— 2nT\x l — y 1 \]. This 
is the expected behaviour at large temperature. 



7 Discussion and Conclusions 

Returning to (5.6) and noting that 

\rj(ir)\ 2 r)(ir)r)(i / V) 



P " w 



L/0) (7.1) 



is left invariant under an exchange of /3 with L, and that the same is true for K(x, x), we 
may view (5.10) as dependence of the zero-temperature chiral condensate on the spatial 
extension L. Thus Figure 1 can be intepreted as the change of {^f^)(L) due to finite 
(spatial) size effects. In particular, if the compactified spatial size shrinks to zero, the 
chiral condensate behaves as 

(**)(£) — >~e-* /Lm i for L^O (7.2) 

and vanishes exponentially (and non-analytically) . 

In the decomposition (5.4) of the bosonic path integral we have integrated over dif- 
ferent L r (l)-principal bundles labelled by the index k. We have assumed that the relative 
normalization between functional integrals with different k is one. To recover the 8- 
vacuum structure [6] we must allow for a relative phase exp(i6) between the sectors k and 
fc + 1 (similarly as in the path integral for a particle on a punctured plane). Then the 
expectation values (%?P±ty) would pick up the phases exp(±i#) and thus the right hand 
sides of (5.10) and (5.12-14) would be multiplied by cos(6). The same happens if one adds 
the explicit ^-dependent term iQ^jl-n to the classical action (1.1). 
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In the conventional path integral solution of the Schwinger model on the plane the zero- 
temperature result (5.13) is derived indirectly, via clustering of the 4-point function (1.17). 
This seems to be necessary since a (naive) direct calculation of the chiral condensate yields 
the wrong result (\t , \[ / ) = 0. This may not come as a surprise since one integrates only over 
potentials with vanishing fluxes and on the other hand we have seen that on the torus only 
configurations with fluxes ±2tc contribute to the expectation value (5.13). Integrating over 
all gauge potentials should yield the correct result (5.13). More precisely, on the plane a 
gauge potential can be decomposed as 



Afj, = A^ - e^dvcf), (7.3a) 



where A is a vortex field 



1 $ 

2^ a 2 + r< 



K = " — -JT~^^ v x\ $ = / E (7.36) 



and SA^ = —e^ v d v (j) a local fluctuation with vanishing flux about the vortex potential. 
One can show that p possesses exactly \k\ zero modes, where \k\ is the integer part of 
|$|/27r [22] (besides these zero modes there is also a resonance state). Since only potentials 
with |/c| = 1 contribute to the chiral condensate one obtains 

/ d$ p($) / V<j> e" r M tr (^P+ip) 

where we have allowed for a Jacobian p of the transformation (5.2) (without harmonics 
h). Again the effective action V contains two pieces, namely a local part coming from the 
chiral anomaly and a global vortex determinant. Since on the plane there aren't different 
topological sectors it would be interesting to see how the ^-parameter emerges in such a 
calculation. 
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